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f>f-\ ' Abstract The secular evolution process, which slowly transforms the morphology of a 

QQ , galaxy over its lifetime, could naturally account for observed properties of the 

l/^ ■ great majority of physical galaxies if both stellar and gaseous accretion pro- 

^^3 ' cesses are taken into account. As an emerging paradigm for galaxy evolution, its 

^3 , dynamical foundation had been established in the past few years, and its obser- 

\l ' vational consequences are yet to be fully explored. The secular evolution picture 

^^ ' provides a coherent framework for understanding the extraordinary regularity 

ri . and the systematic variation of galaxy properties along the Hubble sequence. 
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The dominant view over the past few decades has been that the structural 

properties of galaxies remain largely unchanged unless galaxies were perturbed 

by violent events such as mergers (Toomre & Toomre 1972). This view is 

C^ ■ particularly favored by the currently popular hierarchical clustering/cold dark 

matter (CDM) paradigm of structure formation and evolution (Peebles 1993 
and the references therein). In the early 1980s, photometric and kinematic ev- 
idence in the buldges of late-type galaxies, and hints from N-body simulations 
of barred galaxies which incorporated a dissipative gas component, prompted 
several investigators to speculate that a fraction of late-type bulges might be 
formed from gas accretion under the bar potential (Kormendy 1982; Combes 
& Sander 1981). These early observations and speculations have since been 
further developed into one version of the secular evolution scenario, which 
emphasizes the role of dissipative gas accretion in the formation of the so- 
called "pseudo bulges" in late type galaxies (Kormendy & Kennicutt 2004 and 
the references therein). 

Even though the role of mergers during the early phases of galaxy assembly 
is still to be assessed, growing evidence has shown that at least since z ~ 1 



the rate of merger appears to have been significantly reduced (Conselice et al. 
2003), and during subsequent time the significance of merger in transforming 
galaxy morphology is likely to be overshadowed by the slower secular evolu- 
tion process (Kormendy & Kennicutt 2004). 

However, secular evolution involving gas alone under barred potential leads 
to some apparent paradoxes. First of all, as emphasized by Andredakis and 
coworkers (Andredakis, Peletier, & Balcells 1995), the continuity of the galaxy 
properties across the entire Hubble sequence, highlighted for example by the 
continued variation of Sersic index n in fitting the bulge surface density pro- 
file, indicates that there is not an apparent break in the formation mechanism 
between the late-type disk galaxies and mid-to-early-type disk galaxies. How- 
ever, due to the paucity of gas compared to stars in most galaxies, manifest- 
ing as a gaseous-to-stellar mass ratio of 1/10 or less (which is true even for 
high-redshift late-type disk galaxies, presumably a result of rapid star forma- 
tion following the accretion shocks of dissipative disk formation, D. Sanders 
2002 private communication; F. Combes 2004 private communication), there 
simply is not enough of a reservoir of gas for use to build up the bulges of 
intermediate Hubble type galaxies such as our own, which has a bulge mass 
comparable to the disk mass, not to say for galaxies of even earlier Hubble 
types. Secondly, bulges of intermediate to early type galaxies, including our 
own, consists mostly of stars of very old age (Jablonka, Gorgas & Goudfrooij 
2002), and could not have been built up by the secular accretion of gas over 
a Hubble time which subsequently formed stars. Even though Kormendy and 
Kennicutt (2004) argued that some of these bulges may be termed "classical", 
i.e., they could have formed out of the dissipative collapse at an early stage of 
galaxy formation, they nonetheless pointed out that a good fraction of these 
apparently old bulges have stellar kinematics which were rotation-dominated, 
and are related to the kinematics of their disks, hinting at a secular evolution 
origin for their formation. 

These and other apparent paradoxes, such as the existence in some galax- 
ies of dense core in the central region which appears to be formed by radial 
gas accretion and yet is kinematically decoupled from the bulge stars, can be 
naturally resolved if we allow also the possibility of stellar accretion. Such 
a possibility was not explored in the past decades mainly due to a theoretical 
barrier, i.e., the well-known result that the stellar motion in a galaxy contain- 
ing quasi-stationary non-axisymmetric patterns conserves the Jacobi integral 
in the rotating frame of the pattern (Binney & Tremaine 1987), and the orbital 
motion of stars under such potentials generally will not exhibit secular delay 
or increase, since there is no wave and disk star interaction except at the wave- 
particle resonances which are usually localized for quasi-steady wave patterns 
(Lynden-Bell & Kalnajs 1972). 
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2. Physical Mechanisms for Producing the Secular 
Morphological Evolution of Galaxies 

It was first demonstrated in Zhang (1996, 1998, 1999) that secular orbital 
changes of stars across the entire galaxy disk are in fact possible due to the col- 
lective instabilities induced by the unstable density wave modes such spirals, 
bars, as well as by the skewed 3D density distributions reflected in the twisted 
isophotes of many high-redshift galaxies. These skewed global patterns were 
shown to lead to a secular energy and angular momentum exchange process be- 
tween the disk matter and the wave pattern, mediated by a local gravitational 
instability, or a collisionless shock, at the potential minimum of the pattern 
(Zhang 1996). The integral manifestation of this process is an azimuthal phase 
shift between the potential and density spirals, which results in a secular torque 
action between the wave pattern and the underlying disk matter. As a result of 
the torquing of the wave on the disk matter, the matter inside the corotation 
radius loses angular momentum to the wave secularly and sinks inward. The 
wave carries the angular momentum it receives from the inner disk matter to 
the outer disk and deposits it there, causing the matter in the outer disk to drift 
further out. 

Closed form evolution rates can be derived for this process. In the quasi 
steady state, the rate of angular momentum exchange between a skewed den- 
sity wave pattern (spiral, bar, etc.) and the basic state of the disk, per unit area, 
is given by (Zhang 1996, 1998) 

which, for two sinusoidal waveforms, is given by 

-_(r) = {m/2)A^Avsin{m(t)o), (2) 

at 

where A^, and Ay are the amplitudes of the density and potential waves, re- 
spectively, and 00 is the phase shift between these two waveforms. 

The orbital decay rate it induced on a disk star (or a gas clump) can be 
derived from equation (2) as 

— = --F2wotan(i)sin(m(/)o), (3) 

where F is the fractional wave amplitude (which is the geometric average of 
the fractional density and potential wave amplitude), vq is the circular velocity 
of the star, i is the pitch angle of the spiral, m is the number of spiral arms, 
respectively. 

From the above expressions we can easily see why the skewness of the pat- 
tern is needed: It is the skewness of the mass distribution which leads to the 



potential/density phaseshift through the Poisson equation (Zhang 1996). With- 
out phaseshift, there will be no secular angular momentum exchange between 
the wave and the disk matter at the quasi-steady state, and thus no secular mass 
redistribution. A perfect oval or triaxial mass distribution without any skew- 
ness will not be the direct driver for secular evolution. 

For our own Galaxy, if we assume that over the past Hubble time the Galac- 
tic spiral pattern has an average 20% amplitude and 20 degree pitch angle, the 
orbital delay rate can be calculated using equation (3) to be about 2 kpc of or- 
bital decay per Hubble time. Therefore, a star in the Sun's orbit will not make it 
all the way in to the inner Galaxy in a Hubble time. However, the correspond- 
ing mass accretion rate across any Galactic radius inside corotation is about 
6 X 10^ Mq per Hubble time. A substantial fraction of the Galactic Bulge can 
thus be built up in a Hubble time. 

Another important consequence of spiral-induced wave-basic state interac- 
tion is the secular heating of the disk stars. Since a spiral density wave can only 
gain energy and angular momentum in proportion to the pattern speed i7p, and 
a disk star loses its orbital energy and angular momentum in proportion to its 
circular speed Q, an average star cannot lose the orbital energy entirely to the 
wave, the excess energy serves to heat the star. The resulting orbital velocity 
diffusion rate is given by 

D = {n- %)F'^vI tan(i) sin(m0o). (4) 

Using parameters appropriate for the Galaxy, this can quantitatively explain 
the observed age-velocity dispersion relation of the solar neighborhood stars 
(Zhang 1999; Wielen 1977). This secular heating process, since it originates 
from a local gravitational instability at the arms, serves to increase all three 
space velocities of stars. The vertical component of the velocity dispersion al- 
lows the stars to gradually drift out of the galactic plane as they spiral inward 
to become bulge stars. The corresponding energy injection into the interstel- 
lar medium, if used as the top-level energy input for turbulence cascade, can 
quantitatively account for the size-linewidth relation of the Galactic molecular 
clouds (Zhang et al. 2001; Zhang 2002). 

The secular morphological evolution process leads to the Hubble type of an 
average galaxy to evolve from the late to the early (Zhang 1999). The secular 
evolution speed is expected to be faster for cluster galaxies than for isolated 
field galaxies, which were observed to evolve slowly, because the rate of sec- 
ular evolution is proportional to the wave amplitude squared and pattern pitch 
angle squared (see, e.g., equations 3, 4), and cluster galaxies are found to have 
large amplitude and open spiral patterns excited through the tidal interactions 
with neighboring galaxies and with the cluster potential as a whole. This natu- 
rally accounts for the morphological Butcher-Oemler effect. 
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The secular evolution of galaxies through the internal process mediated by 
skewed non-axisymmetric patterns is a rigorous consequence of the Newtonian 
dynamics and the assumption of global self-consistency. The analytical equa- 
tions for the evolution rates (equations 3 and 4) are quantitatively confirmed 
in N-body simulations (Zhang 1998). Detailed comparison of the predictions 
of the secular evolution theory with the observed galaxy properties and with 
other theories of galaxy evolution can be found in (Zhang 2003). In what fol- 
lows, we will highlight the role of secular evolution in generating the structural 
properties and scaling relations of galaxies. 

3. Generation of the Structural Properties of Galaxies 

Along the Hubble Sequence through Secular Evolution 

The properties of disk galaxies vary systematically along the Hubble se- 
quence from the late to the early Hubble types. An important descriptor of 
the structural characteristics of disk galaxies is the rotation curve. After as- 
sembling hundreds of observed rotation curves for nearby disk galaxies, Per- 
sic, Salucci, & Stel (1996) found that these rotation curves fall onto a two- 
dimensional surface in the three-parameter space of normalized galactic radius, 
velocity, and the absolute magnitude. Since for the nearby galaxies the varia- 
tion of the absolute magnitude corresponds to the variation of Hubble type (i.e., 
bigger and brighter galaxies generally have earlier Hubble types), this means 
that the rotation curve shapes also vary systematically for galaxies of varying 
Hubble types. In Figure 1, we plot three typical rotation curves generated using 
the data from Persic et al. (1996). 

Using dynamical relations for equilibrium axisymmetric/spherically sym- 
metric mass distributions, i.e. V{R)'^ = '^'^ , as well as M^yniR) = 

/q Yjproj{r)2Tirdr = ^q /9(r)47rr^dr, where Md2;„(i?) is the dynamical or to- 
tal mass within a galactic radius R (including the contributions from both the 
luminous and the dark components), and T,proj{R) is the total projected sur- 
face density (disk plus spheroidal) at radius R, we can infer the underlying 
mass and surface density distributions corresponding to the different regimes 
for each type of rotation curves (c.f. Figure 1). 

(I) Steeply rising solid body rotation curve, as observed for the inner regions 
of early type galaxies. 

This regime is characterized by solid body rotation ((7 = JIq = constant, or a 
nearly linear rise of rotation velocity V{R) = QqR). Using V{R)'^ = ^qR^ = 
GMdyn/R, we obtain Mdyn = ^^}^ , which imphes a constant volume density 
P — Po = -Q^- Therefore the projected surface density can be obtained by 
equating Mdyn{R) = j^^proj{r)2'Krdr and Mdyn{R) = |vri?^po, which 
gives T,proj{R) = "^PqR, i-e. the solid-body regime corresponds to a constant 
column density po and linearly rising projected surface density "EprojiR) oc R. 




Figure 1. Typical rotation curves for disk galaxies of varying absolute magnitudes, as well 
as four characteristic regimes which are analyzed in the main text. Top: M=18; middle: M=21; 
bottom: M=24, where M is the absolute magnitude of the representative galaxy (data from 
Persic etal [1996]). 



(II) Gently rising rotation curve, as representative of the very late type disks, 
including most of the dwarf galaxies and low-surface-brightness (LSB) disk 
galaxies. 

This regime can be shown to correspond to approximately constant pro- 
jected surface density, as found to be the case for many LSB disks. Assume 
^proj{R) = So = constant, we obtain Mdyn{R) = Jq T.o2'Krdr = Y.qttR'^, 
which we equate to V{RfR/G to obtain ^(i?)^ = ^/T^ottRG, which indeed 
corresponds to a gently rising rotation curve. 

Instead of constant projected surface density, many LSBs are found to have 
instead constant volume density (de Blok 2003). These two views (constant 
projected surface density, or constant mass volume density) can be consistent 
if the scale height in the central regions of these LSBs are nearly constant, as 
for a pancake-type proto galaxy cloud collapse. 

(III) Flat rotation curve, as is representative of the outer part of most Sb/Sc 
galaxies. 

In this regime, the projected surface density profile should be close to 1/R, 
as shown below. Assume T,proj{R) = "j^" , where Sq = ^proj{R = -Ro)> 
we have M^yniR) = /cf T.proj{'f)27rrdr = T,qRq2ttR. Using again V{Ry = 
GMdyn/R, we obtain ViRf = GT.qRq2txR/R, or V{R) = VGT,oRo2tt = 
constant. Therefore the flat rotation curve regime corresponds to 1/R projected 
surface density distribution. 

(IV) The falling rotation curve regime, as representative of the outer part of 
the early type galaxies. 

In this regime the projected surface density T,proj{R) falls off faster than 
I/R, and can be up to T,proj = (Keplerian). For the extreme case of Keplerian 
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rotation curve due to a concentrated central mass distribution, "Eproj (Router) ~ 

0, and V{Rf = ^fM^^ or V{R) = \f^^^, which is indeed falUng. 

In the secular evolution scenario a typical galaxy starts its life with proper- 
ties similar to that of an LSB galaxy. The high-z correspondence of the local 
LSBs appears to be the damped L^ system (Wolfe 2001), which also possesses 
disk-like morphology with large scale length, and a constituting mass compo- 
nent which is gas rich and metal poor. 

As the proto galaxy disk contracts and condenses, the global spiral insta- 
bility develops, which begins the mass-redistribution process that transforms 
the flat surface density distribution to a more and more centrally concentrated 
mass distribution. The intermediate stage of this evolution resembles Free- 
man's type II disks (Freeman 1970), and the later stage with a a 1/R projected 
surface density resembles a Freeman's type I disk. 

As the evolution progresses, the inward-accreted disk stars rise out of the 
disk plane to become bulge stars. The thermalized three-dimensional bulge 
mass will have nearly constant volume density and undergoes solid-body rota- 
tion. The outer part of the visible disk starts to be drained of matter, producing 
the falling rotation curve. 

This entire evolution sequence therefore corresponds to a galaxy evoling 
from the lower rotation curve to the upper rotation curve in Figure 1 . 

The twisted isophotes often found for the Hubble Deep Field galaxies are 
expected to evolve into the observed variation of major axis with radius for 
elliptical galaxies (Faber & Jackson 1976), possibly bypassing the disk for- 
mation phase, at least for some galaxies. Like the spiral structure, the twisted 
isophotes are results of dissipation and differential rotation, and just as in the 
case of spiral structure, these three-dimensional skewed structures will also 
lead to a phaseshift between the potential and density distributions, resulting 
in accelerated secular morphological changes. 

4. Generation and Evolution of the Galaxy Scaling 
Relations Along the Hubble Sequence 

Galaxy scaling relations (Faber & Jackson 1976; TuUy & Fisher 1977; Djor- 
govski & Davis 1987; Dressier et al. 1987) own their existence first and fore- 
most to the fact that galaxies are equilibrium configurations and satisfy the 
Virial theorem relation (see, e.g., Binney & Tremaine 1987). This relation 
dictates that during any quasi-equilibrium evolution process of a gravitational 
system, —E ~ T ~ —V/2 where T is the kinetic energy, V is the potential 
energy, and E = T + V i& the total energy of the system. The fact that spiral 
galaxies can remain on the TuUy-Fisher relation during secular morphological 
evolution despite the fact that stars can't radiate about 1/2 of the converted po- 
tential energy away, is due to the fact that spiral density wave actually carries 
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about 1/2 of the dissipated orbital energy (equal to potential energy) away to 
the outer disk to be deposited there (Zhang 1999). The existence of a collective 
dissipation process at the spiral arms also made possible the regulated conver- 
sion of orbital velocity to random velocity, so a single kinetic energy (the half 
that is not being carried away) can enter the Virial relation no matter how it is 
apportioned between the circular and the random motions. 

According to the Virial theorem (we set G=l in the following derivations 
of scaling relations) Mdyn = V^Re, where Mdyn is the dynamical mass of 
a galaxy, Ve is the effective velocity spread which corresponds to the total 
velocity dispersion for an elliptical galaxy or the maximum rotation velocity 
for a spiral galaxy, and Re is an effective radius which makes the Virial relation 
hold. 

Define the mean surface brightness of a galaxy SB as SB = -Mj where L 
is the luminosity of the galaxy, and the dynamical mass-to-light ratio M^yn/L 
can be written as M^yn = {^dynl L) * L. Therefore it follows that 

(Mdyn/Lf * l2 = Miy^ = V^R^ = y/^, (5) 

or 

\2 - ^' 



{Mdyn/LYL = ^. (6) 



Therefore we finally have from 

^ = ^'^SBiM,yjL)^- ^^^ 

We thus see that in order to obtain the classical Tully-Fisher relation, which has 
L oc V^, we must have SB ■ (Mdyn/L) r^ constant. This can be accomplished 
in two ways: to have SB and M^yn/L each being constant; or, to have the 
variations of the two factors offset each other during an evolution process. 

Traditionally, the former was assumed to be the case (see, e.g. the discus- 
sions and references in Shu 1982). It is now known that both these quantities 
vary considerably for galaxies along the Hubble sequence. The surface bright- 
ness is found to be higher for earlier type galaxies (McGaugh & de Blok 1998 
and the references therein), whereas the dynamical mass-to-light is found to be 
lower for the earlier Hubble types (Zwaan et al. 1995; Bell & de Jong 2001). 
The observed opposing trends of variation of surface brightness and dynamical 
mass-to-light ratio is naturally explained as the outcome of the secular baryonic 
mass accretion and the increase fraction of luminous baryon mass, especially 
towards the central region of a galaxy as its Huuble type evolve from late to 
early. The Tully-Fisher relation can be maintained as long as the increase in 
surface brightness is accompanied by a corresponding decrease in the dynam- 
ical mass to Ught ratio during the secular evolution process. 
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In the secular evolution scenario, since most of the elliptical galaxies (i.e. 
those lower luminosity, so-called disky ellipticals) are the end results of evo- 
lution from the initial condition of disks, we speculate that spiral galaxies may 
satisfy similar kind of fundamental plane relation just as ellipticals. This has 
been found to be so (see, e.g., Pharasyn, Simien, & Heraudeau 1997). The 
fundamental plane relation for spirals can likewise be derived from the Virial 
theorem. Starting from 

K' = G^, (8) 



we have 



L Rp Ri L 



which leads to 

10 log V; = -(1 + 2/3)Mi - ^le + constant, (10) 

where Mj is the absolute magnitude and ^g is the face-on average surface 
brightness in mag/arcsec^, and where we have assumed yidyj^ oc L'^. Pha- 
rasyn et al. (1997) found that fitting I band and K band data of a sample of 
spiral galaxies to this relation resulted in /3 « —0.15, i.e. the mass-to-light ra- 
tio slightly decreases with increasing luminosity, which is consistent with what 
we would expect from a secular evolution picture, i.e. as the secular evolution 
advances, L increases and Mrfj/„/L decreases. 

The fundamental plane relation for elliptical galaxies and bulges were first 
obtained by Faber et al. (1987), Djorgovski & Davis (1987), and Dressier et 
al. (1987). One type of the fundamental plane relation for ellipticals can be 
written as (c.f. equation 2a and 2b of Djorgovski & Davis 1987) 

L ~ v;3-455i?-°-86. (11) 

This relation can be re-written into the form 

Mt{Re) = -8.62(log Ve + O.l^e) + Constant, (12) 

which can be further written as we arrive at 

10 log K = -1.25M(i?e) - yWe + Constant. (13) 



Comparing equations (10) and (13), we see that apart from a possible differ- 
ence in the "constant" term (accounting for the change in total luminosity), the 
only difference between the spiral and elliptical fundamental plane relations 
is in the different sign of the mass-to-light ratio exponent (3: (3 = —0.15 for 
spirals and (3 = 0.13 for ellipticals and bulges, where j3 is defined through 
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Mdyn/L = L^. This difference is apparently brought about by the fact that 
spiral galaxies still have varying reserves of baryonic dark matter to form stars, 
therefore as the secular evolution proceeds (and therefore L increases) the 
mass-to-light ratio decrease; whereas elliptical galaxies have essentially ex- 
hausted their central baryonic dark matter supply, thus the ellipticals in more 
advanced stage of evolution (which also generally have larger L) will experi- 
ence greater degree of dimming, which is reflected in the increase of mass-to- 
light ratio with L. 

This research was supported in part by funding from the Office of Naval 
Research. 
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